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ABSTRACT 



We present a study of the surface density profiles of the clusters of galaxies 
from the Palomar Distant Cluster Survey (Postman et al. 1996). The survey 
contains a total of 79 clusters of galaxies, covering the estimated redshift range 
of 0.2 ^ z ^ 1.2. We have analyzed the richest clusters in this sample and 
find that the typical Palomar cluster has a surface density profile of r^^'^ 
(r > 0.10 Mpc) and a core radius of 0.05 Mpc. There may be an 
indication that the slope of the surface density profile steepens with increasing 
redshift, though the observational uncertainty is at present too large to be 
conclusive. Our cluster population is inconsistent at a 99.9% confidence level 
with a population of azimuthally symmetric clusters. 

Subject headings: galaxies: clustering; cosmology: observations 



1. Introduction 



Clusters of galaxies provide a powerful probe of the nature of galaxy formation and the 
origin of large scale structure in the universe. Having a well-studied sample of intermediate 
and high-redshift clusters is vital to understanding the evolution of galaxies and large 
scale structure. Because the diversity of cluster properties and the effects of evolution are 
significant at high redshift (Gunn & Dressier 1988; Bower et al. 1994; Castander et al. 
1994; Postman et al. 1996), it is essential to sample many distant clusters with different 
properties. Large scale simulations have shown that cluster properties such as profile shape 
and substructure can constrain the nature of large scale structure formation theories and 
the mass density of the universe (Evrard et al. 1993; Crone, Evrard & Richstone 1995; Jing 
et al. 1995; Tyson & Fischer 1995). Detailed studies of these observational parameters have 
already been made in some individual, intermediate redshift clusters. Using weak lensing 
to map the mass distribution in several clusters. Small et al. (1995) have studied two X-ray 
luminous clusters 1455+22 {z = 0.26) and 0016+16 {z = 0.55), Tyson & Fischer (1995) 
have studied Abell 1689 {z = 0.18), and Squires et al. (1995) have studied Abell 2218 
{z — 0.175). All four clusters show moderate to extreme degrees of structure. However, 
the distribution of structure in mass, light, and X-rays are all well correlated at radii 
larger than 100 kpc. The projected mass density (and light) profile of A1689 can be 
well approximated by a power-law of while the projected total mass, gas, and 
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light surface densities of A2218 are consistent with an isothermal sphere (r~^), though 
perhaps shghtly steeper. Smail et al. (1995) found that, while the galaxies arc very good 
tracers of the mass, they are less concentrated, the respective core radii for 1455+22 being 
^mass _ 50+25 ^pc and r^"' = 90^25 kpc. A multi-color photometric study of three 
clusters, A3284 {z = 0.15), A3305 {z = 0.16), and A1942 {z = 0.23), also yield consistent 
core radii of r^'*' ~ 120, 100, and 120 h^^ kpc, respectively (Molinari et al. 1994). 

In this paper, we expand on these detailed analyses of individual intermediate redshift 
clusters by studying a large sample of intermediate and high redshift clusters from the 
Palomar Distant Cluster Survey (hereafter PDCS; Postman et al. 1996). We examine both 
the individual and global properties of the cluster profiles, specifically the cluster profile 
slope, core radius, and degree of asymmetry. Because the PDCS uses a completely objective 
and automated algorithm to detect clusters, it provides us with the largest, statistically 
complete sample of distant clusters. In addition, the selection biases due to the detection 
technique can be well quantified through simple Monte-Carlo simulations. 

In §2 of this paper we briefly describe the PDCS and the cluster sample used in this 
analysis. The parameters of the individual cluster profiles are discussed in §3. Composite 
cluster profiles as a function of redshift are presented in §4. Cluster morphology is discussed 
in §5. We present a preliminary comparison of our results to large-scale cosmological 
simulations in §6 and summarize the results in §7. 



2. The Cluster Sample 

The cluster sample is derived from an optical/near IR survey with the 4-shooter CCD 
camera on the Palomar 5 meter telescope. The cluster catalog is the subject of the first 
paper in this series (Postman et al. 1996; hereafter Paper I). However, we discuss briefly 
the aspects of this catalog which are necessary for the following analysis. 

The Palomar Distant Cluster Survey was conducted in two broad band filters, the 
F555W and F785LP of HST's Wide Field/Planetary Camera. We denote these bands V4 
(F555W) and I4 (F785LP) according to the convention in Paper I. The area covered consists 
of five one-square degree areas around the sky, all at galactic latitudes of 30 degrees or 
greater. The data are complete to = 23.8 and I4 — 22.5. A matched filter algorithm was 
used to objectively identify the cluster candidates by using positional and photometric data 
simultaneously. This technique is likely to be more robust than previous optical selections 
which simply looked for surface density enhancements, a method which can be significantly 
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affected by superposition effects {e.g., Gunn, Hoessel & Oke 1986; Coudi et al. 1991). An 
advantage of tliis teclinique is tliat redsfiift estimates of the cluster candidates are produced 
as a byproduct of the matched fiher; the main disadvantage is that we must assume a 
particular form for the cluster luminosity function (for the flux filter) and cluster radial 
profile (for the radial filter). The radial filter P(r) and the fiux filter L{m) are given by 



P(r) = , ^ - , ^ if r < rco 

otherwise 



(hfm—m*) in-0-4(m-m*) 
^ 0[m m) iU 

b{m) 

P{r) is an azimuthally symmetric cluster surface density profile which has a characteristic 
core radius (rc) and which falls off at large radii as r~^. The function is explicitly cut off at 
an arbitrary cutoff radius (rco). We have chosen Vc — 100 kpc and Tco = 10 x Tc such 
that the radial profile resembles the profiles of nearby clusters and that we have optimized 
the cluster detections relative to the spurious detection rates (Paper I). (f){m — m*) is the 
differential Schechtcr luminosity function with a = —1.1 and M* = —21.0 and —21.9 in the 
V4 and bands, respectively; b{m) is the background galaxy counts. For the derivation of 
these filters and a detailed explanation, see Paper I. We have used extensive simulations in 
both Paper I and this paper to quantify the selection bias due to the functional form of 
the matched filter. We find that this selection bias has a minimal effect on the resulting 
distribution of profiles. The matched filter algorithm does a good job at preserving the true 
distribution of profile shapes for a broad range of cluster profile parameters. 

Candidate clusters are detected individually in each band and then matched with each 
other to locate those systems which are detected in both bands. The catalog consists of 
79 clusters of galaxies detected with estimated redshifts between 2; ~ 0.2 and 1.2; the 
uncertainty in the estimated redshift {zest) is cr^^^t ~ 0.2 (see Paper I; Lubin 1995). 87% 
of the cluster candidates are matched detections (detected in both the V4 and I4 bands). 
The amplitude of the matched filter provides an estimate of the cluster richness. This filter 
richness (A^./) is a measure of the effective number of L* galaxies in the cluster (see §4.2.2 
of Paper I; Lubin 1995). Through Monte-Carlo simulations, we can statistically determine 
the relation between Ac; and the actual cluster richness as determined by the specification 
of Abell (1958). This relation is dependent on profile shape; as the cluster profile slope 
steepens, a given Ac; value corresponds to a lower richness class. For a cluster which has 
a surface density profile of approximately (the assumed profile of the radial filter). 
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h-ci ^ 40 corresponds to Abell R > 1 (Paper I). Therefore, in the analysis which follows, we 
will examine only those clusters with > 40. The cluster sample that will be the subject 
of this paper consists of 57 clusters in the V4 band and 50 clusters in the /4 band. Table 1 
lists the ID number of these clusters (according to the convention of Paper I). 



3. Individual Cluster Profiles 



3.1. The Profile Fittings 



We begin by examining the individual profiles in the V4 and 1^ bands. We have 

an estimate of the cluster redshift {zest) and the cluster center from the matched filter 
detection algorithm. The cluster center is taken to be the position of the peak cluster 
signal for each detection (see Paper 1). Surface density profiles of each cluster candidate 
are computed by counting the number of galaxies in successive annuli around the 
center of each cluster. We bin out to a given physical radius of 1.0 h^^ Mpc (assuming 
Q,o — ^iHo — 100 h km s~^ Mpc~^, and h — 0.75) at the estimated redshift of the cluster. 
Twenty bins of equal radial extent are used. The uncertainty in each bin is calculated 
from Poisson statistics. An annulus of radii 1.0 < r < 1.5 Mpc around each cluster is 
used to measure directly the background level for each cluster. The measured background 
varies ~ 10 — 15% from cluster to cluster and field to field. As an example, at a redshift of 
z ~ 0.5, the number of cluster members (out to a radius of 1.0 h^^ Mpc) in a R = 1 cluster 
is ~ 150 galaxies or 10 - 20% of the total number of galaxies in that region. 

We quantify the slope at large radii (described by a power-law exponent a) and the 
cluster core radius (rc) of the cluster profile by fitting each resulting background-subtracted 
surface density profile to a King model (King 1962; Cavaliere & Fusco-Femiano 1976; 
Sarazin & Bahcall 1977) : 



S{r) 



1 + 



-3/3+1 



(3) 



where So and Tc are the normalization and core radius of the King model, respectively. (3 
parameterizes the density fall-off at large radii; that is, the surface density profile falls off 



as r 



-3^+1 



Fig. 1 shows the distribution of power-law exponent {a) and core radii (rc) for those 
cluster profiles which had acceptable fits {'X^/v < 1.5; u — 17). The average la single 
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parameter errors (Avni 1976) in the fit parameters of an individual cluster profile are 
Aa/a and Avc/r^ =-60%^° i^-d-' Table 3). We are able to constrain the slope 

of the surface density profile much more accurately than the cluster core radius (see also 
§3.2). As a consistency check, we have examined the effect on the profile fittings by (1) 
calculating the cluster background at a larger limiting radius; (2) changing the bin size (by 
±50%); and (3) using the weighted center rather than the peak center of each cluster (where 
the weighted center is defined as the mean position as weighted by the amphtude of the 
matched filter signal within the cluster detection; see also §5); the cluster fit parameters are 
all consistent within the measurement error with Aa/a ^ 10% and Atc/tc ^ 50%. 

The variation in the cluster profile parameters is large, as is the sensitivity to the 
background level. Therefore, it is difficult to say anything quantitative on an individual 
basis (see also §4). Table 2 lists the median values and "Icr" confidence interval for the 
distribution of a and Tc in each band. We define a la confidence interval around each 
parameter by determining the points in the distribution such that the area in each tail is 
16%. The average PDCS cluster has a surface density slope of (a) ^1.4 and core radius 
of (tc) ~ 0.05 h^^ Mpc. Note that the distribution of core radii in the Palomar clusters 
indicates that approximately 20 - 30% of the clusters have very small or essentially no core 
radii (Fig. 1). The median values in the V4 and I4 bands are consistent with each other 
(Table 2). 

The observed range in a is completely consistent with a sample of nearby, rich clusters 
{{a) ^ 1.6; Bahcall 1977; see also Beers & Tonry 1986; Oegerle et al. 1985; Postman, Geller 
& Huchra 1988), the galaxy-cluster cross correlation function which represents the average 
net galaxy density proffie around rich clusters ((a) a; 1.2 — 1.4; Seldner & Peebles 1977; 
Peebles 1980; Lilje & Efstathiou 1988), and measures of the mass (and light) distribution of 
intermediate redshift clusters (Small et al. 1995; Tyson & Fischer 1995; Squires et al. 1995). 
The range of core radii of the PDCS clusters is also typical of other well-studied nearby 
and intermediate redshift clusters. A study of 15 rich, nearby clusters yields core radii of 
Tc = 0.12 ± 0.02 Mpc (Bahcall 1975; though a similar study of 12 clusters by Dressier 
1978 indicate core radii larger by a factor of ~ 2). Small et al. (1995) examined the core 
radii of two intermediate redshift. X-ray luminous clusters and found Tc = 0.09lo;o3 h^^ Mpc 
for 1455+22 {z = 0.26) and Vc ^ 0.16 h'^ Mpc for 0016+16 {z = 0.55), while Molinari et al. 
(1994) found Tc ^ 0.10 - 0.12 h'^ Mpc for A3284 {z = 0.15), A3305 {z = 0.16), and A1942 
(^ = 0.23). 

In Fig. 2, we show the cluster profiles for three matched PDCS cluster candidates 
(cluster ID # 003, 036 and 059; Paper I). In order to show the diversity of the PDCS 
clusters, we have chosen clusters which cover a range in estimated redshifts (zest) and profile 
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shape. Table 3 lists the relevant parameters of the fittings to these cluster profiles. 



3.2. Evaluating Selection Biases 



In order to estimate the measurement error and to demonstrate that the matched filter 
does not significantly bias the profile distribution, we have performed simple Monte-Carlo 
simulations. We estimate the errors in our fitting procedure and evaluate the selection 
function by running the cluster finding algorithm on artificial fields which match, as closely 
as possible, the galaxy and cluster distributions of the actual PDCS fields. We reproduce 
our matching criteria in the cluster detection method by creating ten pairs of galaxy fields, 
each 1 deg^, as viewed through the V4 and I4 passbands. We accomphshed this by first 
modehng the actual background density - magnitude function of regions of the PDCS 
fields which contain no detected clusters. For each simulated field, we randomly distribute 
the appropriate number of background galaxies in each magnitude bin according to this 
function. On top of this uniform field, wc superimpose 25 simulated clusters with the 
condition that the separation between clusters must be greater than 500 arcscc. The clusters 
are randomly given redshifts between 0.2 and 1.2 (in z increments of 0.1) and richnesses of 
R > 0. The only restriction placed on cluster richness is that clusters at z < 0.5 cannot be 
richer than richness class 2 and clusters at 2; > 0.6 cannot be poorer than richness class 
1. We have chosen for the cluster luminosity function the same Schechter function as the 
matched filter; however, we assume that half of the clusters contain spirals (and apply a 
Scd spiral k-correction) and half contain ellipticals (and apply an elliptical k-corrcction). 
The k-correction will affect the redshift estimate of the cluster (cr^^^^ ^ 0.2) but not the 
richness estimate Ad (Lubin 1995; Paper I). 

Firstly, we estimate the measurement error inherent in our fitting procedure by 
simulating a population of clusters with a unique surface density slope (a) and core radius 
(rc). We choose the average profile parameters observed in the PDCS : a = 1.4 and 
Tc = 0.05 h^^ Mpc. The cluster galaxies arc distributed in an azimuthally symmetric 
manner according to the appropriate King profile (Eq. 3). We apply the same detection 
criteria and richness cut (Ac; > 40) to the simulated clusters as for the PDCS clusters. For 
each simulated cluster, a surface density profile is created and then fit to the King model in 
the manner described in §3.1. The resulting distributions of fit parameters for 20 simulated 
clusters are shown in Fig. 3. The median values and their la errors (as defined in §3.1) 
of these distributions are listed in Table 4. The median values of the distributions of fit 
parameters reproduce well the actual input parameters of the simulated clusters; however. 



-8- 



the distributions are artifically broadened by ~ 10% in a and ~ 50% in r^- These errors 
are due to uncertainties, most notably background subtraction, associated with the fitting 
procedure; in addition, errors in the estimated redshift contribute to the broad width of 
the Tc distribution. The measurement uncertainties associated with this particular cluster 
population are typical of the errors for cluster populations with other combinations of a 
and Tc- 

Secondly, we examine the selection bias of the radial filter [a = 1; Tc = 0.1 hr"^ Mpc) by 
exploring a broad range in cluster shape through two profile distributions. In the first set of 
simulations, we generate a Gaussian distribution of profile shapes which is reasonably close 
to that observed in the PDCS (Table 2); the simulated clusters have a surface density slope 
distribution of a = 1.4 ± 0.6 and a core radius distribution of — 0.05 ± 0.05 hT^ Mpc. In 
the second set of simulations, we generate a Gaussian distribution of profile shapes which 
is very different from the PDGS, a = 2.5 ± 0.5 and = 0.10 ± 0.05 hr^ Mpc. As above, 
the cluster galaxies are distributed in an azimuthally symmetric manner according to the 
appropriate King profile (Eq. 3); the same detection criteria and richness cut (Ac/ > 40) 
are applied. Fig. 4 shows the input parameter distributions (shaded histograms) of the 
simulated clusters and the output parameter distributions (solid line histograms) of those 
clusters detected by the algorithm. There appears to be no significant selection bias in 
either of the two simulated cluster populations due to the matched filter; we reproduce 
the input distributions of a and over this broad range of cluster parameters extremely 
well. As noted in Paper I, clusters with steeper profiles yield the strongest matched filter 
signal. This will affect the cluster detection probability, though the exact effect is also 
dependent on the cluster core radius and the cluster richness. However, we have shown with 
the simulations of this section that, for reasonable choices of the input cluster distribution 
(that is, distributions with standard deviations typical of that observed in the PDCS), the 
matched filter does not significantly bias the detected (output) parameter distribution. 

We have quantified both the measurement error associated with the profile fitting 
procedure, as well as the selection bias associated with the functional form of the matched 
filter. The observed median parameters of the PDCS profile distributions appear to provide 
an accurate indication of the true parent population of cluster shapes. 



3.3. Central Density — Richness Comparison 

We can also estimate the slope of the cluster surface density profiles by examining the 
relation between richness {N^) and central density {N^ of the individual clusters. Similarly 
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to Abell (1958), central density and richness are defined as the number of member galaxies 
(above background) that are brighter than + 2™ (where mo, is the magnitude of the third 
brightest galaxy) and located within a projected radius of 0.25 Mpc and 1.0 Mpc, 
respectively (see §4.2.3 of Paper 1 for details). Fig. 5 shows Nr versus No for the PDCS 
clusters in the two passbands (top panels). The median values of Nr/No are 2.5l];;o (^4 
band) and 2.7^1'^ {I4 band). The standard deviation in these distributions represent la 
confidence intervals as defined in §3.1. 

The expected relation between richness and central density for the simulated cluster 
population described in §3.2 (a = 1.4 ± 0.6 and = 0.05 ± 0.05 h^^ Mpc) is shown 
in the bottom panels of Fig. 5. The median values of Nr/N^ are 2.4lJ { (V4 band) and 
2.7li;o (-^4 band), consistent with those of the PDCS. We note that the relation between 
Nr and No for a population of clusters with a much steeper surface density profile {e.g., 
r~^; {Nr/Ng) ~ 1.6) is inconsistent with that observed in the PDCS; the probability, as 
measured by the KS test, that the PDCS distribution is drawn from the same parent 
population as the Nr/Ng distribution of an cluster population is ^ 10~^^%. 



4. Cluster Profile Evolution 

In an effort to improve our constraints on the cluster profiles, we explore the average 
profile shape, specifically the surface density slope a, as a function of redshift by creating 
global composites of our sample of PDCS clusters. Because we exclusively use estimated 
redshifts (observed redshifts are presently available for fewer than 10 PDCS clusters) whose 
uncertainty are estimated at a^^^^ ^ 0.2 (see §4.2.1 of Paper I), we sort each cluster by 
its estimated redshift (zest) in Va and I4 bands into three broad redshift ranges : (1) 
0.2 < Zest < 0.4, (2) 0.5 < Zest < 0.7, and (3) 0.8 < Zest < 1-2. We examine the cluster 
profile at r > 0.10 Mpc by creating a composite cluster profile. We logarithmically bin 
galaxies in the radial range 0.1 < r < 1.0 Mpc around each cluster center with a bin 
width of 0.1 in log r (i.e. 10 bins). The background level determined for each individual 
cluster profile (from an outer annulus; see §3.1) is first subtracted. We stack the cluster 
profiles in each redshift interval according to their physical radius (0.1 < r < 1.0 Mpc). 
The resulting composite cluster profiles are then fit to a power-law of the form S{r) — Sgr""'. 



4.1. Composite Cluster Profiles 
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The composite profiles of the PDCS in the V4 and I4 bands are shown in Fig. 6. The 
redshift range is indicated in the upper right-hand corner of each paneL The error in the 
mean value is determined from the variance in each radius bin. The solid line indicates the 
best-fit power-law function. Table 5 lists the number of cluster profiles used to make the 
composite, the best-fit power-law slopes a, and their la single parameter errors. These 
errors are calculated according to the procedure described by Avni (1976) who presents 
the correct "minimum x^" method for calculating confidence limits in a fitting with one 
"interesting parameter" (i.e. a). The slopes of the composite profiles are consistent with the 
results of the individual fittings (Table 6) . The surface density slopes of the most distant 
clusters {zest > 0.8) are shghtly steeper than their intermediate redshift counterparts. This 
may simply reflect the larger uncertainty in the fitted parameters when co-adding such a 
small number of clusters (9 in V4 band and 5 in I4 band in this redshift interval). Such a 
large dispersion is also seen in the simulations when co-adding a similar number of high 
redshift clusters. In addition, it may be the result of a bias due to the inherent difficulty 
in subtracting an appropriate background for clusters at very large redshift. Nonetheless, 
some evolution in the cluster profile is expected for certain cosmologies (see §6). 

Since richer clusters are believed to have formed at earlier epochs (Peebles 1993), we 
may expect to see a variation in profile shape with richness. We, therefore, examine the 
composite cluster profiles of a richer subset {Aa > 70) of the PDCS clusters. We find that, 
within the observational uncertainty, there is no difference between the power-law slopes of 
the composites of this richer population and the composite profiles presented in Fig. 6. 



4.2. Evaluating Selection Biases 

In order to demonstrate that the composite profiles are not being biased by a few 
clusters, we performed a parametric bootstrap analysis (see Kendall & Stuart 1967). This 
technique involves randomly drawing a new sample (of equal number) from the original 
cluster sample. From each new sample, we create a composite cluster profile as described 
above and fit the power-law function. This procedure was performed on 50 resampled 
datasets. The average (and standard deviation) of the resulting a values are completely 
consistent (within la) with the fits to the original true composite in all three redshift 
intervals (Table 5). 

We confirm that our technique for creating cluster composites is not biased due 
to either the method of background subtraction or the functional form of the matched 
filter by performing two consistency checks. Firstly, we make composite cluster profiles 
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from regions of the PDCS fields which contain no detected clusters. For each "blank" 
region, wc randomly assign an estimated redshift drawn from the actual estimated redshift 
distribution. We then bin the galaxies, subtract the background, and stack the profiles in 
the same manner used to create the PDCS composites. The resulting composites are shown 
in Fig. 7. All profiles are consistent with zero, indicating that our method of background 
subtraction is reasonably accurate. 

Secondly, we use the simulations described in §3.2 to examine whether our cluster 
composites accurately refiect the mean chistcr profile as a function of redshift. In an 
identical manner to the PDCS, we have created cluster composites from those clusters 
detected in 10 simulated fields. We have examined the simulated cluster population which 
closely resembles the PDCS clusters {a = 1.4 ± 0.6 and Tc = 0.05 ± 0.05 h'^ Mpc). The 
results are shown in Fig. 8. The best-fit a parameters are hsted in Table 6. The resulting 
composite profiles of the simulated clusters are consistent with the mean slope in each 
redshift interval and in each band. 

We have shown that our actual cluster composites (Fig. 6) are a reliable representation 
of the global cluster population and are not significantly biased, as a function of redshift, 
by the intrinsic nature of the matched filter. 



5. The Cluster Morphologies 

Observational studies of nearby clusters have shown that a significant portion of 
clusters exhibit evidence of ellipticity, multiple components, and other forms of substructure 
{e.g., Geller & Beers 1982; Dressier & Schectman 1988; West & Bothun 1990; Jones & 
Forman 1992; Mohr et al. 1995; Buote & Canizares 1995; de Theije, Katgart & van Kampen 
1995). Presently, measurements such as the axial ratio, centroid shift, and ellipticity have 
been used to quantify the cluster "morphology" (the degree of asymmetry on large scales) 
in the cluster surface density and X-ray contours (Mohr, Fabricant & Geller 1993; Mohr et 
al. 1995; Buote & Canizares 1995; de Theije et al. 1995). 

We take a similar approach by examining the degree of symmetry in the PDCS clusters 
through the "filtered" images of the PDCS fields, i.e. the amplitude of the matched filter 
signal as a function of position within a given cluster detection. At each position, this 
signal measures how accurately a cluster centered at that position matches the filter. Since 
the filter is an axisymmetric, monotonically decreasing function of radius, the amplitude as 
a function of position is effectively a map of the structure in the cluster candidate. The 
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first two columns of Fig. 9 show the contour plots of sample PDCS clusters in each of the 
three redshift intervals (top to bottom). Adjacent panels show the contours of each cluster 
in the V4 and /4 band, respectively. The side of each panel corresponds to 2 h^^ Mpc at the 
estimated redshift of the cluster. The contours represent lines of constant amplitude of the 
matched filter signal. 

We compare these contour maps with the azimuthally symmetric, centrally concentrated 
clusters of our simulations (§3.2). The last two columns of Fig. 9 show the contour levels 
of simulated clusters at the same estimated rcdshifts as the PDCS clusters. The simulated 
symmetric clusters may be marginally rounder than the PDCS clusters. We quantify 
any difference by measuring the offset between the peak center and the weighted center 
of the cluster. The peak center (used in the profile analysis) is defined as the position 
of the maximum signal in a given cluster detection; the weighted center is defined as 
the mean position as weighted by the amplitude of the matched filter signal within the 
cluster detection. The magnitude of this offset is a measure of the degree of asymmetry 
of the clusters. Fig. 10 shows the distribution of the offset (in kpc) for the PDCS 
clusters (shaded histograms) and for the azimuthally symmetric simulated clusters (solid 
line histograms). Using the KS test, the null hypothesis that the offset distributions are 
drawn from the same parent cluster population can be rejected at a 99.9% confidence level 
for the V4 band and at a greater than 99.9999% confidence level for the /4 band, indicating 
that the typical PDCS cluster is not azimuthally symmetric. That is, the PDCS clusters are 
inconsistent with an azimuthally symmetric cluster population with its characteristic a and 
core radius Tc- Significant structure on large scales, including ellipticity and bi-modality, 
in both mass, light, and X-rays have been previously found by Small et al. (1995), Tyson 
& Fischer (1995) and Squires et al. (1995) in the intermediate redshift clusters 1455+22 
{z = 0.26), 0016+16 {z = 0.55), AbeU 1689 {z = 0.18), and Abell 2218 {z = 0.175). 



6. Compcirison to Structure Formation Simulations 



Recently, Crone et al. (1994) and Jing et al. (1995) have used high resolution N-body 
simulations to examine the dependence of cosmology on the cluster mass density profile 
and morphology. The PDCS provides a large statistical sample to which we can compare 
the results of these cosmological simulations. Crone et al. (1994) examine the dependence 
of the cluster profile on both cosmology and initial density field by using a scale-free initial 
power spectra P{k) oc k"' with initial spectral indices n = —2, —1 and for four cosmologies 
: Einstein deSitter {fig — 1), open {flo — 0.2 and 0.1), and flat, low-density {flo — 0.2, 
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Ao = 0.8). Jing et al. (1995) examined the Standard CDM model and six low-density CDM 
models with and without a cosmological constant. The resulting average mass density 
profiles are all well-fit by a power-law p{r) oc for radii greater than ~ 0.2 h^^ Mpc (or 
local density contrasts between 100 and 3000). There exists a clear trend toward steeper 
slopes with decreasing Qg- Crone et al. (1994) also found a correlation between steepening 
slopes and increasing spectral index n in their initial power spectrum. In addition, Jing 
et al. (1995) find that, for a given ^o{< 1), clusters in the flat model (Aq = 1 — Op) have 
flatter mass density proflles than in the corresponding open model. A composite of these 
results are presented in Fig. 11 where we plot the slope of the surface mass density proflle 
{a — 7 — 1) versus for simulated clusters at 2; = 0. Flat cosmologies are indicated by 
closed circles, while open cosmologies are indicated by open circles. 

We would now hke to compare the composite surface density proflles of the PDCS 
clusters with the cluster simulations; however, the simulations model the total mass 
distribution, whereas we examine the galaxy distribution. It is now possible to map the 
total mass density distribution as a function of radius through observations in X-rays and 
through gravitational lensing of rich clusters of galaxies {e.g., Henry, Briel & Nulsen 1994; 
Small et al. 1995; Tyson & Fischer 1995; Squires et al. 1995). These observations reveal 
that there is excellent agreement between the projected structure of the galaxy, gas, and 
dark matter distributions on scales larger than ~ 100 kpc. Therefore, a comparison 
between the composite surface density proflles of the PDCS clusters with those expected 
from the large scale simulations should allow reasonable constraints to be placed on the 
mass proflles. As discussed in Tyson & Fischer (1995), our mean projected density slope 
of {a) 1.4 is intermediate between standard CDM and an Qo = 0.35 open CDM or A+ 
CDM simulated profiles but appears to be a better match to the latter two. In Fig. 11, 
we indicate by two dotted lines the range of power-law slopes that are obtained from the 
composite PDCS cluster profiles over the rcdshift range 0.2 ^ z ^ 1.2 (Table 5 and Fig. 6). 
The flo = I profiles are indeed on the low end of our observed a distribution; however, this 
is far from conclusive. 

Up to this point there has been little examination of the evolution of cluster profile 
shape as a function of redshift because the cluster samples at intermediate to high redshift 
have been sparse. The Palomar survey provides the first statistically complete cluster 
survey over a large redshift range. The cluster composites may indicate that the proflles get 
steeper with increasing redshift. Such a trend in the mean slope is most prominent in the I 4 
band though the trend is not signiflcant given the large observational uncertainty. However, 
such an effect may also be observed in simulated clusters from high-resolution (20 kpc) 
N body simulations (Xu 1995). Because clusters form much later in an flo — ^ universe, 
standard CDM shows the largest evolution in the cluster proflle. The slope decreases from 
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a = 1.51 ± 0.27 = 1) to q; = 1.29 ± 0.19 (-z = 0). The shallowing of the power-law slope 
at radii greater than ~ 0.2 hr^ Mpc is largely the result of an increase in the cluster core 
radius due to continued merging between z = 1 and z = (Xu 1995; Pen 1995). Profiles of 
clusters in an open {Qo = 0.35) CDM shows no detectable evolution with a = 1.38 ± 0.38 
{z = 1) and a = 1.38 ± 0.19 (2; = 0). + A models yield intermediate results. Other 
non-standard cosmologies, such as mixed dark matter models, also show a steepening of 
the cluster profile with redshift (Walter & Klypin 1995). At present, with httle redshift 
information, it is not possible to use the observed variation to rule out any cosmological 
model. 

The cluster morphology, or the degree of symmetry, has also been used to delineate 
between cosmological models (Evrard et al. 1993; Mohr et al. 1995; Jing et al. 1995). 
N-body simulations and hydrodynamic simulations show that the surface density and X-ray 
image contours of clusters in low-density CDM models {e.g., Qg — 0.2) are much more 
regular and centrally concentrated that those in an Einstein-de Sitter model [fig — 1). 
Clusters in low cosmologies are older and are, therefore, more relaxed and symmetric 
than clusters in an = 1 universe. Mohr et al. (1995) have used these results to argue in 
favor of a fiat universe, given the observed degree of asymmetry and structure in present-day 
clusters. Jing et al. (1995) found, however, that low- flo models with a cosmological constant 
A {e.g., Qo = 0.3 and A = 0.7) also produce a large fraction of clusters with significant 
deviations from spherical symmetry, though the density and X-ray image contours are still 
rounder than the Qg — ^ case. However, this matter is still controversial as Pen (1995) 
and Xu (1995) find little quantitative difference in the structure characteristics of clusters 
formed in these cosmological models. We have used the contours of constant amplitude of 
the matched filter signal to examine the degree of symmetry in the PDCS clusters. We do 
find that the PDCS clusters are inconsistent with a population of azimuthally symmetric 
simulated clusters (Fig. 10). 



7. Summciry 

We have examined the galaxy distributions of the richest clusters of galaxies from the 
Palomar Distant Cluster Survey. Through an analysis of individual cluster profiles and 
composite profiles as a function of redshift, we find that the typical Palomar cluster has a 
profile of r~^'^ at radn greater than 0.10 Mpc and a core radius of 0.05 Mpc. The 
distribution of core radii in the Palomar clusters indicates that up to 30% of the clusters 
have very small or essentially no core radii. Using simple Monte-Carlo simulations, we have 
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shown that the median cluster profile parameters appear to be an accurate representation of 
the actual cluster population; the measurement errors associated with these parameters are 
~ 10% in the surface density slope a and ~ 50% in the cluster core radius r^. The average 
profile parameters of the PDCS clusters are consistent with measures of nearby clusters, 
as well as other intermediate redshift clusters (Molinari et al. 1994; Tyson & Fischer 1995; 
Squires et al. 1995; Small et al. 1995). There may be an indication that the mean cluster 
profile steepens with increasing redshift, though the observational uncertainty is presently 
too large to be conclusive. In addition, we find that a large fraction of the PDCS clusters 
have a significant degree of asymmetry. The PDCS cluster population is inconsistent with 
a population of circularly symmetric clusters of galaxies at a greater than 99.9% confidence 
level. 
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Figure Captions 

Figure 1 : Distribution of surface density slopes (a) and core radii (rc in h^^ Mpc) from 
the fittings to the PDCS cluster profiles. The shaded and solid line histograms 
represent the distribution of parameters from individual profile fittings of the PDCS 
clusters for the V4 and /4 bands, respectively. We show only those fittings which have 
acceptable x^/^ (see §3.1). 

Figure 2 : King profile fittings to three matched cluster candidates. Adjacent panels 
show the cluster profile in the V4 and /4 bands, respectively. The estimated redshift 
{zest) and the cluster ID number is indicated in the top right and left of each panel, 
respectively. A solid line indicates the best-fit King model (Eq. 3). 

Figure 3 : Distribution of surface density slopes (a) and core radii (rc in Mpc) from 
the fittings to the profiles of a simulated cluster population with a discreet a and Tc 
combination (a = 1.4; Tc = 0.05 Mpc). The median values accurately reproduce 
the actual profile parameters, though the distributions are artificially broadened due 
to inherent errors in the fitting process (see §3.2). 

Figure 4 : The input (shaded histograms) and output (sohd line histograms) distributions 
of surface density slopes (a; top panels) and core radii (rc in h'^ Mpc; bottom panels) 
of the two sets of simulated clusters (see §3.2). The input parameter distributions are 
specified in the upper right-hand corner of each panel. All distributions are shown 
separately for the V4 and I 4 bands. 

Figure 5 : Abell richness Nr (< 1.0 Mpc) versus central density No (< 0.25 Mpc) 
for the PDCS clusters and for simulated clusters (see §3.2) in the V4 and I4 bands. The 
black dots, open squares, and stars indicate clusters detected at estimated rcdshifts of 
0.2 < Zest < 0.4, 0.5 < Zest < 0.7, and Zest > 0.8, respectively. The expected relation 
from the simulations is consistent with the observations. 

Figure 6 : Composite surface density profiles (0.10 < r < 1.0 Mpc) of the PDCS 
clusters. The clusters have been divided in three redshift bins based upon their 
estimated redshifts in each band. The individual background surface density has been 
subtracted from each individual cluster profile before creating the composite. The 
variance in each bin is used to compute the error in the mean value. The best-fit 
power-law functions are shown. 

Figure 7 : Composite surface density profiles (0.1 < r < 1.0 h"^ Mpc) of regions in 
the PDCS fields which contain no detected clusters. We have assumed estimated 
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redshifts such that we have roughly the same number of profiles in each redshift bin 
as the real cluster composites (Fig. 6). The constant background surface density 
has been subtracted from each "profile" before creating the composite. All profiles 
are consistent with zero, indicating that our method of background subtraction is 
reasonably accurate. 

Figure 8 : Composite surface density profiles (0.1 < r < 1.0 h'^ Mpc) of the simulated 

clusters. The clusters have been divided in three redshift bins based upon their 
estimated rcdshifts in each band. The constant background surface density has been 
subtracted from each individual cluster profile before creating the composite. The 
variance in each bin is used to compute the error in the mean value. The best-fit 
power-law functions are shown. 

Figure 9 : Contour maps of the amplitude of the matched filter signal in both bands 
for one cluster in each of the three redshift intervals (top panel 0.2 < z^st < 0.4; 
middle panel 0.5 < z^st < 0.7; bottom panel 0.8 < Zest < 1-2). The first two columns 
show actual PDCS clusters, and the last two columns show azimuthally symmetric, 
simulated clusters. The side of each panel corresponds to 2 Mpc at the estimated 
redshift of the cluster. 

Figure 10 : Distribution of offsets (in h^^ kpc) between the peak center (the position 
of the maximum signal of the matched filter) and the weighted center. The shaded 
histograms represent the distribution of offsets of the PDCS clusters. The solid 
line histograms represent the distribution of offsets of the azimuthally symmetric, 
simulated clusters. The distributions are shown separately for the V4 and I4 bands. 
The actual data and the simulations are statistically different (see §5). 

Figure 11 : The expected power-law slopes of the cluster mass surface density profiles 
versus flo from various cosmological simulations (Crone et al. 1994; Jing et al. 1995). 
Open and closed circles indicate open and flat cosmologies, respectively, and represent 
simulated clusters at 2; = 0. Dashed lines indicate the range of power-law slopes of the 
composite PDCS cluster galaxy profiles over our full redshift interval 0.2 ^ 2; ^ 1.2 
(Table 5). 
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Table 1 : The Sample of PDCS Clusters 



Band PDCS Cluster ID Number 

001 002 003 004 006 008 010 Oil 012 015 
016 017 018 019 020 021 022 023 024 025 

027 030 031 033 034 035 036 037 038 039 
041 042 044 047 048 049 050 051 052 053 
054 055 056 057 059 062 063 065 066 067 
068 069 071 072 075 076 078 
001 002 003 004 006 008 009 010 Oil 012 
014 015 016 017 018 019 020 021 022 023 
024 031 033 034 035 036 037 039 041 042 
045 049 050 051 052 054 055 056 057 059 
061 062 063 064 067 068 069 075 076 079 



Table 2 : Median Values of Surface Density Slope a and Core Radius Tc {h ^ Mpc) of the 

PDCS Clusters 



Band 


# 


a 




h 


43 
37 


i.OO_o.47 
±.00_o.35 


o.o5l°:S^ 
o.o5l°:J° 



Table 3 : Surface Density Slope a and Core Radius Tc {h ^ Mpc) of Three Matched PDCS 

Clusters 





PDCS 




h 




Panel 


ID# 


a Tc 


a 




Top 


036 


1 79+0.39 n r\Q+OM 
-■-•'^-0.23 U.Ui)_o.02 


9 co+0.69 
^•00-0.42 


0.21+°1? 


Middle 


003 


1 7Q+0.66 n qc+0.05 


1 90+0. 18 
-'-•^'-'-0.13 




Bottom 


059 


^■-■-•J-O.SQ U.U4_o.o2 


9 49+0.51 
^•^^-0.30 


0.041°:°^ 
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Table 4 : Measurement Error Associated with Median Values of Surface Density Slope a 
and Core Radius Tc {h~^ Mpc) from the Simulated Clusters 



Band 


a 




h 


1.0U_Q 23 

1 4Q+0.12 
-■-■^•5-0.15 


0.0571°:°?! 
0.059t°:°P 



Table 5 : Power-Law Slopes of Composite Surface Density Profiles of the PDCS Clusters 



Redshift Interval 




^4 




/4 




# 


a 


# 


a 


0.2 < Zest < 0.4 


26 


1 cin+o-22 

-'-•3U_o.34 


21 


1 00+0.13 
-L-'J^-O.IS 


0.5 < Zest < 0.7 


22 


1.34l°i? 


24 




0.8 < Zest < 1-2 


9 


2.06l°i^ 


5 


l-93l°i! 



Table 6 : Power-Law Slopes of Composite Surface Density Profiles of the Simulated 

Clusters 



Redshift Interval 








h 




# 


a 


# 


a 


0.2 < Zest < 0.4 


29 


J--^0_Q 20 


29 


1 47+0.15 

-0.19 


0.5 < Zest < 0.7 


66 


-, /iQ+0.15 


71 


J-'^^-o.ie 


0.8 < Zest < 1-2 


45 


1 07+0.25 
-'-•'^'-0.55 


53 


1 44+0-15 

-■-•^^-0.20 
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Fig. 1. — Distribution of surface density slopes (a) and core radii (rc in Mpc) from the 
fittings to the PDCS cluster profiles. The shaded and solid line histograms represent the 
distribution of parameters from individual profile fittings of the PDCS clusters for the V4 
and I4 bands, respectively. We show only those fittings which have acceptable x^/^ (see 
§3.1). 
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Fig. 2. — King profile fittings to three matched cluster candidates. Adjacent panels show 
the cluster profile in the V4 and I4 bands, respectively. The estimated redshift (-Ze^t) and the 
cluster ID number is indicated in the top right and left of each panel, respectively. A solid 
line indicates the best-fit King model (Eq. 3). 
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Fig. 3. — Distribution of surface density slopes (a) and core radii (r^ in Mpc) from the 
fittings to the profiles of a simulated cluster population with a discreet a and Vc combination 
(a = 1.4; Tc = 0.05 Mpc). The median values accurately reproduce the actual profile 
parameters, though the distributions are artificially broadened due to inherent errors in the 
fitting process (see §3.2). 
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Fig. 4. — The input (shaded histograms) and output (sohd hne histograms) distributions of 
surface density slopes {a; top panels) and core radii {r^ in Mpc; bottom panels) of the 
two sets of simulated clusters (see §3.2). The input parameter distributions are specified in 
the upper right-hand corner of each panel. All distributions are shown separately for the V4 
and I 4 bands. 
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Fig. 5. — Abell richness Nr (< 1.0 Mpc) versus central density No (< 0.25 Mpc) 
for the PDCS clusters and for simulated clusters (see §3.2) in the V4 and I4 bands. The 
black dots, open squares, and stars indicate clusters detected at estimated redshifts of 
0.2 < Zest ^ 0.4, 0.5 < Zest < 0.7, and Zgst > 0.8, respectively. The expected relation 
from the simulations is consistent with the observations. 
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Fig. 6. — Composite surface density profiles (0.10 < r < 1.0 Mpc) of the PDCS clusters. 
The clusters have been divided in three redshift bins based upon their estimated redshifts 
in each band. The individual background surface density has been subtracted from each 
individual cluster profile before creating the composite. The variance in each bin is used to 
compute the error in the mean value. The best-fit power-law functions are shown. 
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Fig. 7. — Composite surface density profiles (0.1 < r < 1.0 Mpc) of regions in the 
PDCS fields which contain no detected clusters. We have assumed estimated redshifts such 
that we have roughly the same number of profiles in each redshift bin as the real cluster 
composites (Fig. 6). The constant background surface density has been subtracted from 
each "profile" before creating the composite. All profiles are consistent with zero, indicating 
that our method of background subtraction is reasonably accurate. 
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Fig. 8. — Composite surface density profiles (0.1 < r < 1.0 Mpc) of the simulated 
clusters. The clusters have been divided in three redshift bins based upon their estimated 
redshifts in each band. The constant background surface density has been subtracted from 
each individual cluster profile before creating the composite. The variance in each bin is 
used to compute the error in the mean value. The best-fit power-law functions are shown. 
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Fig. 9. — Contour maps of the amplitude of the matched filter signal in both bands for 
one cluster in each of the three redshift intervals (top panel 0.2 < Zest < 0.4; middle panel 
0.5 < Zest < 0.7; bottom panel 0.8 < Zest < 1-2). The first two columns show actual PDCS 
clusters, and the last two columns show azimuthally symmetric, simulated clusters. The side 
of each panel corresponds to 2 Mpc at the estimated redshift of the cluster. 
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Fig. 10. — Distribution of offsets (in kpc) between the peak center (the position of 
the maximum signal of the matched fiher) and the weighted center. The shaded histograms 
represent the distribution of offsets of the PDCS clusters. The solid line histograms represent 
the distribution of offsets of the azimuthally symmetric, simulated clusters. The distributions 
are shown separately for the V4 and I4 bands. The actual data and the simulations are 
statistically different (see §5). 
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Fig. 11. — The expected power-law slopes of the cluster mass surface density profiles versus 
Qo from various cosmological simulations (Crone et al. 1994; Jing et al. 1995). Open 
and closed circles indicate open and flat cosmologies, respectively, and represent simulated 
clusters at z = 0. Dashed lines indicate the range of power-law slopes of the composite 
PDCS cluster galaxy profiles over our full redshift interval 0.2 ^ z ^ 1.2 (Table 5). 



